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ABSTRACT 

Motivated by recent observations of short-timescale variations in the infrared 
emission of circumstellar disks, we propose that coronal mass ejections can remove 
dust grains on timescales as short as a few days. Continuous monitoring of stellar 
activity, coupled to infrared observations, can place meaningful constraints on the 
proposed mechanism. 

1. Introduction 

Two recent papers show that the radiation associated with warm, dusty circumstellar 
disks around cool stars can vary on surprisingly short timescales. 

Melis et al. (2012) report on observations of the star TYC 8241 2652 1 taken at around 
10 iim and 20 //m in which the flux dropped by a factor of about 30 between 2008 May 6 
and 2010 Jan 8-10, a period of less than two years. TYC 8241 2652 1 (hereafter TYC 8241) 
is thought to be a star of approximate spectral type K2, and an age ~10 Myr, situated at a 
distance of about 140 pc. 

Meng et al. (2010) report on 24 /im variability in two solar-like stars ID8 and HD 23514. 
In ID8 the flux was observed to increase by about 30 percent over a two year period, and to 
then decrease by a similar amount over about a year. In HD 23514 the flux was observed 
to drop by about 10 percent over a period of a few months, or less. ID8 is a G6-G7 dwarf 
at a distance of ~450 pc, and HD 23514 is of spectral type F5V at a distance of ~130 pc 
(e.g., Soderblom et al. 1993). Age estimates for both stars are in the range 30-130 Myr. 

Motivated by these observations, we examine in the current paper potential mechanisms 
for obtaining short-timescale variability. In particular, we propose that coronal mass 
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ejections (CMEs) may play a hitherto unrecognized role in producing variability in the IR 
excess of debris disks. 

In Section 2 we describe the observations and some of the suggested interpretations, 
and we discuss CMEs and their potential effects. In Section 3 we briefly discuss dust grain 
properties, and our conclusions follow. 



2. The IR Fluxes and Variability 

Mehs et al. (2012) model the pre-2009 IR excess with emission from an optically 
thin collection of grains of size 0.3 //m and a temperature of 450 K, orbiting TYC 8241 
at a distance of RiO.4 AU. They find that the ratio of IR excess to stellar luminosity is 
Lm/L^ 0.11. They also note that given the luminosity of TYC 8241 (L* ^ 0.7 Lq), grains 
with radii of ~0.2 /im or smaller would be radiatively ejected from the system. Based on 
these parameters, they estimate that roughly 5 x 10^^ g in grains with radii of order 0.3 /im 
are required to produce the pre-2009 epoch IR excess. 

Similar values apply to the stars reported by Meng et al. (2010), except that in ID8 
and HD 23514 the fractional IR luminosities are lower (Lir/L* ^ 0.02), and the stars are 
more luminous, implying that the dust has to be placed somewhat farther away, at around 
~1 AU. Meng et al. (2010) also note that the spectra of HD 23514 and ID8 (and of a few 
other extreme excess stars) indicate the presence of fine crystalline silicates or silica dust, 
requiring grain sizes of order 0.1 /xm. This is well below the blowout size (see below) for 
these stars. 

The grain blowout size due to radiation pressure is given by (see e.g., Chen & Jura 
2001) 

^ 167rM,pc 
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where L* and M* are the stellar luminosity and mass, Q is the radiation pressure coupling 
coefficient, p is the density of an individual grain, and c is the speed of light. Grains whose 
radii satisfy condition (1) are expelled by radiation. For grains much larger than the mean 
photon wavelength (as expected for the stars under discussion), the effective cross section 
can be approximated by the geometrical cross section, so that Q f« 1. Consequently, for 
grain density p ~ 2 g cm~^ we find for solar parameters that ~ 0.3 iim.. 

For the types of stars we are considering here, the assumed grain sizes are close to the 
blowout limit. Consequently, we can assume that the IR emission is associated with the 
smallest grains that can survive being blown out. 

2.1. Conditions for Increasing the Flux 

Only one of the three stars that exhibits short-timescale mid-IR variabihty, IDS, shows 
an unambiguous increase in the flux (by about 30% over a period of less than two years; 
Meng et al. 2010). This increase, however, may not require the occurrence of any particular 
event, but rather it may simply represent the replenishment of the small-scale end of a 
standard cascade distribution in the sizes of particles (of the type {n{D) oc D~^-^). In 
other words, if the disk were to start with the entire coUisional cascade in place (with 
a distribution of sizes from dust to large bodies), then the IR flux would decrease if a 
mechanism were found to remove dust grains of sizes ~1 /im or smaller. Following such a 
decrease, however, the timescale to regenerate the observationally inferred mass of ~0.1 ^m 
particles by collisional cascade is only of order 0.1 yr — much shorter than the observed 
timescale for the flux increase (Meng et al. 2010; Melis et al. 2012). 

Consequently, we will concentrate in the present work on identifying a mechanism that 
can decrease the flux on a timescale of less than a year, by way of removal of dust at the 
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short end of the size distribution. 

2.2. A Mechanism for Decreasing the Flux 

Mehs et al. (2012) considered the possibihty of vaporizing the small grain population 
using the high-energy flux from a stellar flare. They found, however, that x-rays are 
rather ineffective at heating grains much larger than 0.3 /im. Furthermore, they showed 
that to obtain the necessary vaporization would have required a total x-ray energy of 
E 10^^ ergs, and an x-ray luminosity of 10^-10^ times the stellar luminosity. Mehs et al. 
(2012) noted that flares of such magnitudes have never been observed. 

Here we would hke to examine the potential effects of coronal mass ejections (CMEs), 
which for the Sun, occur together with solar flares in solar eruptive events. These are 
occasions when magnetic reconnection causes a reconfiguration of the magnetic field and 
liberation of energy, which goes into mass motions, particle acceleration and plasma heating. 
First, we note that the observed total kinetic energy released in solar CMEs is comparable 
to the radiated flare energy (Emslie et al. 2005, 2012). This is consistent with theoretical 
expectations from models of magnetic reconnection through magnetohydrodynamic 
processes. The latter suggest that energy is released equally into ohmic heating and into 
fluid motions (e.g.. Priest & Forbes 2000). About 25% of the bolometric flare energy in 
solar flares arises in the form of UV and soft x-ray radiation (Benz 2008). 

We can estimate how large a CME would be required to remove the grains that 
contribute to the IR excess, by taking the CME to consist of an ionized plasma which is 
coupled to a magnetic held. As long as the dust grains are charged and have sufficiently 
small gyroradii (see discussion in Section 3), they would be swept up by such a magnetic 
plasma. The dust grains would be charged by the ionizing radiation of the flare that 
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accompanies the CME. 

The ejection velocity of the mass in the CME, Vcme, can be expected to be of the order 
of the escape velocity from the central star, therefore 

The velocity distribution of solar CMEs is actually quite broad, ranging from a few tens of 
km to 3000 km (Chen et al. 2006). While the velocity could be somewhat reduced 
(e.g. by interplanetary draf; Gopalswamy et al. 2006), there is a wide range of values that 
would give short timescales for dust removal (see equation (4) below). 

If we assume that the mass of dust that needs to be expelled is mdust — 10^^ g, and that 
the dust is located at a radius of i?dust — 1 AU, then from the conservation of momentum 
we need the mass in the CME, mcME, to satisfy 

"^dust 

mcME = rUdustTf — - 0.07mdust ■ (3) 

where V/"*** is the escape velocity of the dust (at i^dust)- Thus, we require mcME ~ 10^° g 
in order to expel a mass of mdust ~ 10^^ g. We note that while CMEs are somewhat 
directional, the distribution of their opening angles extends all the way to about 160 
degrees, with more energetic CMEs having larger opening angles (e.g., Gopalswamy et al. 
2009; Aarnio et al. 2011). While CMEs have not been conclusively detected on stars other 
than the Sun (Leitzinger et al. 2011), the largest solar flares have a nearly 100% association 
with large CMEs (Yashiro et al. 2006). The precise relation between mcME and the flare 
energy is not known, but estimates based on equipartition between flare-radiated energy 
and CME kinetic energy (e.g., Emslie et al. 2005, 2012; Osten et al. in prep.) give (for 
mcME ~ 10^° g) a required flare energy of (a few) x 10^^ ergs, within the realm of previously 
observed large stellar flares. This calculation envisages a single flare/CME, but in principle, 
several smaller flare/CME pairs occurring within a short period of time (< 10^ s) could also 
do the job. 
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The timescale on which the dust can be removed is given approximately by 

Tremoval ^ = 2.4 X 10^ S — — — — , (4) 



which leads, for the values indicated, to timescales of a few days. 



3. Grain Properties 

In order for a CME (an ionized magnetic plasma) to be able to sweep up the dust 
grains, those grains need: (i) to be charged, and (ii) their gyroradii need to be sufficiently 
small, Tgyro < -Rdust) for them to be able to react to the flowing CME. Otherwise, any 
substructure in the CME should not affect the process in any significant way. The first 
condition appears to be easily satisfied. Since the work function for grains is of order 5 eV 
(e.g., Tielens 2010), we need to consider the number of photons with energies > 10 eV, i.e. 
the FUV and XEUV wavelength ranges. In a flare with an energy of -Eflare ~ 10^^ ergs, one 
can expect the number of such photons to be A^ph ~ 6 x 10^^. At the same time, the number 
of grains is of order Ng = rudust/fng ~ 10^*^ (where nig is a typical mass of a grain). Since 
the electron yield per photon is of order Ig/la ~ 0.1 (Tielens 2010), where /„ is a typical 
photon absorption length and le is a typical electron collision length, and Np^ ^Ng, we can 
safely assume that the net charge per grain, Z, satisfies Z ^ 1. 

The gyroradius of a grain is given by 



Ze B -4AU(^J (^2.3xl0-i3<^J U-2xl06cms-V VgJ ' 
where c is the speed of light, e is the electron charge, B is the magnetic field strength (at 
1 AU), and the mass of the grain rrig is scaled with the value for a grain of radius 0.3 /im. 
While magnetic fields have not been measured in stellar flares or CMEs, we can estimate 
the field in the CME by noting that the energy of the accompanying flare is expected to 
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scale with the magnetic energy density in the flaring volume, E oc B'^ (Schrijver 2011). 
Using solar flares as reference, and assuming a dipole field configuration, we would thus 
obtain (B/IQ-'^G) - {En^JlO'^^ erg)i/^ which for E - lO^^ ergs would yield 5 - a few 
G. We note here that Gopalswamy (2006) found solar CMEs at 1 AU to have an average 
value of ~10~^ G, within a factor of ten of the value obtained by the above scahng. 
Remembering that we also expect Z ^ 1, equation (5) therefore indicates that even with 
a potential reduction in the field value due to opening angle considerations or reduction in 
field strength, rgyro < -Rdust, as required for dust removal by CMEs. 

4. Summary and Conclusions 

We have shown that CMEs could remove IR-emitting dust on timescales of days. This 
of course requires the stars in question to be active, and to have significant flares. We note 
that HD 23514, one of the stars that motivated this research, has not been detected in an 
x-ray survey of the Pleiades (Daniel, Linsky, & Gagne 2002). However, one star that has 
been detected (HIT 1122) actually shows lower rotational velocity than HD 23514. Hence, 
it is not implausible that the latter is an active star. 

More importantly, the present work highlights the need for continuous monitoring of 
these and similar stars, to be able to correctly assess their activity level, in conjunction with 
IR flux variations. Furthermore, active M stars with debris disks would provide for ideal 
targets to test the ideas presented in this paper. The currently existing data on transient 
mass loss that is associated with stellar activity (e.g., Leitzinger et al. 2011) is rather 
insufficient for placing meaningful observational constraints. Another promising avenue 
could be to observe the radio signatures of MHD shocks that are expected to be produced 
as CMEs traverse the outer stellar atmospheres at Super- Alfvenic speeds. These would be 
the stellar equivalents of solar "type 11" bursts (e.g., Gopalswamy et al. 2009). 
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